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Abstract: Exact analytical expressions for the AS = 1 coupling Im G^; in terms 
of observable spectral functions are given. This coupling determines the size of the 
A/ = 3/2 contribution to e' . We show analytically how the scheme- dependence and 
scale dependences vanish to all orders in l/N^ and NLO in 05 explicitly both for 
Qi and Qg. Numerical results are derived for both and Qg from the r-data and 
known results on the scalar spectral functions. In particular we study the effect of 
all higher dimension operators. The coefficients of the leading operators in the OPE 
of the needed correlators are derived to NLO m. as- 



Keywords: Kaon Physics, CP violation, Chiral Lagrangians, QCD. 



Contents 



1. Introduction 2 

2. Overview 3 

3. The Qi and Operators 4 

4. Exact Long— Short-Distance Matching at NLO in as 7 

4.1 The Qr contribution 7 

4.1.1 Q7 Long-distance 7 

4.1.2 Q7 Short-distance 7 

4.2 The Qs Contribution 8 

4.2.1 The Disconnected Contribution 9 

4.2.2 The Connected Contribution 10 

4.3 Sum 11 

5. Bag Pcirameters 12 

6. The n^^((5^) Two-Point Function and Integrals over It 14 

7. The Scalar-Pseudo-Scalar Two-Point Function n^V+pp(<5^) 19 

8. Numerical Results for the Matrix-Elements and Bag Parameters 22 

9. Compcirison with earlier results 25 

9.1 Role of Higher Dimensional Operators 26 

9.2 The Large A^c Limit 28 

10. Conclusions 29 

A. Calculation of the Corrections of O(a^) to the Dimension Six Con- 
tribution to TlljiiQ"^) 31 
A.l Renormalization Group Analysis 31 
A. 2 Calculation of the Constants Cf'^ and F^^^ 33 

A. 3 The constants pij 34 

B. Calculation of the Corrections of C(a^) to the Dimension Six Con- 
tribution to ^^ss+pp"^{Q'^) 35 

B. l Renormalization Group Analysis 35 
B.2 Calculation of the Constants Cf^ and F^^^ 36 



1 



B.3 Calculation of the Pij. 



37 



1. Introduction 

The lowest order SU(3) x SU(3) chiral Lagrangian describing [A^l = 1 transitions 
is given by 

^\AS\=i = C Ge tr (As^u^Qu) + CF^ [Gg tr (As^uX) + G'gtr ( A32X+) 

+ Gsyt'^'^'tr (Aiju^) tr {Amu")] + h.c. (1.1) 

with 

C = -l^VudVus* -1.07 • 10-^ GeV-^ (1.2) 
Fo is the chiral hmit value of the pion decay constant /tt = (92.4 ± 0.4) MeV, 



Aij = uXijU^ {kj)ab = SiaSjb (1.3) 



and U = — exp (i\/2$/Fo) is the exponential representation incorporating the 
octet of hght pseudo-scalar mesons in the SU(3) matrix 

The SU(3) X SU(3) tensor f^'^^ can be found in [1] and Q = diag(2/3, -1/3, -1/3) 
is a 3 X 3 matrix which collects the electric charge of the three light quark flavours. 

In the Standard Model in the chiral limit, e' is essentially determined by the 
value of (Im Ge)/G27 and (Im G^)/{Ke Gg) including Final State Interactions to all 
orders. Therefore the knowledge of Im is of primordial importance. The value 
of this parameter is dominated by the electroweak penguin contributions Qi and Qs 
[2]. 

The paper consists out of two parts. In the first part. Sections 2-5, we discuss 
how the X-boson approach takes care of the scheme-dependence in the chiral limit 
independent of the large Nc expansion we used in our previous work. We also show 
precisely how the needed matrix-elements in the chiral limit are related to integrals 
over spectral functions. This clarifies and extends the previous work on this relation 
[3, 4, 5, 6]. Equation (4.27) is our main result, but we also present the expression in 
terms of the usual bag parameters in Section 5. 

In the second part. Sections 6-10, we present numerical results and compare 
our results with those obtained by others and our previous work. Sections 6 and 7 
describe the experimental and theoretical information on both Im ^^r{Q^) and the 
scalar-pseudo-scalar Im ^^ss+pp{Q'^) spectral functions and give the values of the 
various quantities needed. The comparison with earlier results is Section 9. 
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In addition, in the appendices we derive the NLO in 0:5 coefficient of the leading 
order term in the OPE of the needed correlators in the same scheme as used for the 
short-distance weak Hamiltonian. This coefficient was previously only known in a 
different scheme [7]. 

2. Overview 

This section describes the underlying reasoning elaborated in more detail in the 
next two sections. In particular we use a simplified notation here to allow simpler 
intermediate expressions, but we refer to the full equations of the following sections. 

We start from the effective action derived from the Standard Model using short- 
distance renormalization group methods of the form (Eq. (3.1)) 

rsD = E CMQM . (2.1) 

i=7,8 

This effective action can be used directly in lattice calculations but is less easy to 
use in other methods. What we know how to identify are currents and densities. We 
therefore go over to an equivalent scheme using only densities and currents whereby 
we generate (2.1) by the exchange of colourless X-bosons (Eq. (3.5)) 

rx^J2gi{pic)Xl{q'^iq) (2.2) 

i 

where the coupling constants gi can be determined using short-distance calculations 
only. The result is Eqs. (3.6) and (3.7). At this step the scheme-dependence in the 
calculation of the Wilson coefficients Ci is removed but we have now a dependence 
on Mx and the scheme used to calculate with Fx- 

We then need to evaluate the matrix-elements of (2.2). For the case at hand 
this simplifies considerably. In the chiral limit, the relevant matrix element can be 
related to vacuum matrix elements (VEVs). The disconnected contributions are just 
two-quark condensates. The connected ones can be expressed as integrals over two- 
point functions (or correlators) as given in Eq. (3.8), which we evaluate in Euclidean 
space. The two relevant integrals are Eqs. (4.1) and (4.11). 

Both of the integrals are now dealt with in a similar way. We split them into 
two pieces at a scale fi via Eq. (4.2). The two-point function to be integrated over 
is replaced by its spectral representation, which we assume known. 

The long-distance part of the integral can be evaluated and integrals of the 
type (4.4) and (4.21) remain. 

The short-distance part we evaluate in a somewhat more elaborate way which 
allows us to show that the residual dependence on the X-boson mass disappears and 
that the correct behaviour given by the renormalization group is also incorporated. 
To do this, we split the short-distance integral in the part with the lowest dimensional 
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operator, which is of dimension six for both Q-j and Qs, and the remainder, the latter 
is referred to as the contribution from higher-order operators[8]. 

The dimension six part can be evaluated using the known QCD short-distance 
behaviour of the two-point functions at this order. It is vacuum expectation values 
of dimension six operators over for Q-j and over for times a known function 
of as- The vacuum expectation values can be rewritten again as integrals over two- 
point functions and the resulting integrals are precisely those needed to cancel the 
remaining Mx-dcpendence. For the contribution from all higher order operators we 
again perform simply the relevant integrals over the same two-point functions as 
for dimension six and they are the ones needed to match long- and short- distances 
exactly. 

This way we see how our procedure precisely cancels all the scheme- and scale- 
dependence and fully relates the results to known spectral functions. 



3. The Q'j and Operators 

The imaginary part of G^; is dominated by the short-distance electroweak effects 
and can thus be reliably estimated from the purely strong matrix-elements of the 
I AS" I — 1 effective action below the charm quark mass 




with Im Ci = Uilm T the imaginary part of the Wilson coefficients, r = —Xt/Xu and 
Aj = VidV*g. The Hi coefficients are known to two loops [9, 10] and 

Q7 = {salML 2^<i(^P^^1p)r^ (3-2) 

q=u,d,s 

Qs^isa'^Mh ifqi^nl^l^p, (3-3) 

q=u,d,s 

with {QJfiQ) ^fi;^ = ^7/i(l — {+)15)q- Up to 0(q;|), the and Qg operators only mix 
between themselves below the charm quark mass via the strong interaction. 

The QCD anomalous dimension matrix 7(1^) in regularizations like Naive Dimen- 
sional Regularization (NDR) or 't Hooft-Veltman (HV) which do not mix operators 
of different dimension, is defined as^ {i — 7,8) 

^ = - E rii^) QA^) ; 7(^) = E 7^"^«"(^) (3.4) 

'^^ j=7,8 n=l 

^In a cut-off regufarization one has, on the right hand side, an infinite series of higher dimensional 
operators suppressed by powers of the cut-off. Explicit expressions for the matrices 7-''(z^) are in 
App. A. 
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where a(z/) = as{i')/'n. 

At low energies, it is convenient to describe the AS = 1 transitions with an effec- 
tive action V^d which uses hadrons, constituent quarks, or other objects to describe 
the relevant degrees of freedom. A four-dimensional regularization scheme like an Eu- 
clidean cut-off, separating long-distance physics from integrated out short-distance 
physics, is also more practical. In addition, the color singlet Fierzed operator basis 
becomes useful for identifying QCD currents and densities. The whole procedure has 
been explicitly done in [11, 12] and reviewed in [13, 14]. 

At low energies, the effective action (3.1) is therefore replaced by the equivalent 



q=u,d,s 

+ g8{l^c,---) Yl ^q,8Uqd)L + {-2)^eg{sq)R 

q=u,d,s 



(3.5) 



Here all colour sums are performed implicitly inside the brackets. There is also a 
kinetic term for the X-bosons which we take to be all of the same mass for simplicity. 

The couplings gi are determined as functions of the Wilson coefficients by 
taking matrix elements of both sides between quark and gluon external states as 
explained in [11, 12, 13, 14]. We obtain 



\97{l^c)\ 
Ml 



= Im C7(^r) 



l + a(//c) (777 In^ + Ar 



77 



Im Cs{i^r) [a{nc) Aryg] + 0{a{nR) - a{nc)) 



(3.6) 



and 



s(/^c)P 
Ml 



= Im CsilJiR) 
Im C^{|JiR) 



«(Afc) ( 787 In Ar87 



0{a{iiR) - a{iic)) ■ 



(3.7) 



7^^^ is due to the anomalous dimensions of the two-quark color-singlet densities or 
currents. It vanishes for conserved currents. In our case 733^ = —27^), where 7^' is 
the QCD anomalous dimension of the quark mass in the regularization used in (3.5). 
The values of Arij = (r — r)ij have been calculated in [12]. 

The effective action to be used at low-energies is now specified completely. Notice 
that singlet color currents and densities are connected by the exchange of a colourless 
X-boson and therefore are well identified also in the low energy effective theories, 
and the finite terms which appear due to the correct identification of currents and 
densities. 
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The coupling G^; is defined in the chiral hmit so that we can use soft pion 
theorems to calculate the relevant matrix-elements, and relate them to a vacuum- 
matrix-element^. For the contribution of Q-j and Qs, we obtain 



d'px 1 



(27r)4 pj, - M\ 



Where n^'^(p^) is the following two-point function in the chiral limit [3, 4]: 



(3.8) 



p^'p" - g^^p^ 



(3.9) 



In Eq. (3.8) wc used the chiral limit so SU(3) chiral symmetry is exact. L[R)^ — 
{u'-)^d)L{R) or Li^RY = {d'y'^s)L(R), nf;^(p^) vanishes and Ilg'l_^_pp{p'^) is the two-point 
function 

41pp(/) ^ i I d^ye^^-^(0|r[(5 + iP)(«)(y)(5-iP)(«)(0)] |0) (3.10) 



with 



(3.11) 



The 3x3 matrix A^°^ = -\/2 7/-\/3 and the rest are the Gell-Mann matrices normalized 
to tr (a("U(^)) = 26"^ An alternative form for the last term in (3.8) is, 



|^8(/^c,---)p3i 



{2nY Pi - Ml 



n 



SS+PP 



ipl) 



(3.12) 



with 



^55+Pp(p')=^ ld'ye^y^O\T {dd)L{y){ss)R{Q) |0) 
and {qq)^ji) = ?(1 - (+)75)g- 



(3.13) 



^In the real K — ^ tttt case we would need to evahiatc integrals over strong-interaetion five point 
functions, three meson legs and two X-boson legs. For vacuum matrix-elements this reduces to 
integrals over two-point functions, the two X-boson legs. The same is not possible for Gs and G27 
since the corresponding terms are order and have zero vacuum matrix elements. 
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4. Exact Long— Short-Distance Matching at NLO in as 
4.1 The Qr contribution 

In Euclidean space, the term multiplying \g-j\^ in the rhs of (3.8) can be written as 

with = —q^. We split the integration into a short-distance and a long-distance 
part by 

/■oo rix^ POO 

with Mx » H^- In QCD, Hj^jiiQ^) obeys an unsubtracted dispersion relation 

uUQ^)^rM'J-^^. (4.3) 

Jo TT t + Q"^ 

4.1.1 Q7 Long-distance 

Putting (4.3) in (4.1) and performing the integral up to //^ gives 

-1^ r *4 + t) "^hw ^ " (ii) • 

with the use of the Weinberg Sum Rules [15], Eqs. (6.2). 

4.1.2 Q7 Short-distance 

At large Q'^ in the chiral hmit, ^J^jiiQ^) behaves in QCD as [16] 



nL(Q') -EE ^''by (o|o« „,(o)|o>(.) (4.5) 

n=0 i=l V ' ' 



where O2(^^^3)(0) are dimension 2(n + 3) gauge invariant operators. 
oi'\0) = \ L^iO) R,iO) = ^ (srd)^ (0) (d7.«)^ (0) ; 

of (0) = (5 + 2P)(f')(o) (5 - ^p)(°)(o) - (5 + ^p)^'Ho) (5 - ^p)^'Ho) 

= 3(rfd)L(0)(ss)R(0). (4.6) 

The coefficients C^\u,Q'^) are related to the anomalous dimension matrix defined 
in (3.4). This can be used to obtain the NLO in as part of the coefficient with the 
same choice of evanescent operators as in [9, 17] , calculations of the q;| term in other 
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schemes and choices of evanescent operators are in[7]. Our calculation and results 
are in App. A. At the order we work we only need the lowest order [16] 



-787 



(4.7) 



The values of the coefficients of the power corrections are physical quantities and 
can be determined with global duahty FESR^, [18, 19, 20], 



m=0 i=l •^'^^0 * 



n+2 



TT 



(4.8) 



with n > 0. So is the threshold for local duality"^. At leading order in as only n = m 
survive and we can rewrite the short-distance contribution to (4.1) as 
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= a(/i) In -J^i ^ ^ 



M 



X 



(2vr) 



Ml 



+ 
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-1 



C. 



E-E 



2(7t+3) 



1 + 



Ml, 



(0|O«+3)(0)|0) + O(a^) + 



a[iJ,) m -^^i 



Mx J (27r)4 \Ml 



7^ 35,. nr^(5~) - {jifs^ppm - ^fs+ppm) 



+ 



9 



IGtt^ Jo 
where we have used 

4.2 The Qs Contribution 

In Euclidean space, the term multiplying \g^\'^ in the rhs of (3.8), is 



(4.9) 



Q'^ + Ml 



n 



(0-3) 



T ^^SS+PP 



(4.11) 



"^The specific form (4.8) is only true to lowest order in as due to the hi{Q'^) dependence at higher 
orders. 

discussion of the value of the local duality onset is in Section 6. 
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with 

niVp(Q') - nilpp(Q^) = Ti^ss'^piQ') ■ (4.12) 

This two-point function has a disconnected contribution, corresponding to what is 
usually called the factorizable contribution^. We spht off that part explicitly: 



4.2.1 The Disconnected Contribution 

We have included in (3.7) all the 0{as) logs and finite terms that take into ac- 
count passing the four-quark matrix element from the cut-off /ic regulated X-boson 
effective theory to the MS one. Therefore to the order needed 

(0|5W(0)|0)2 = 3(0|gg|0)|l^(^c) (4.14) 



and from now on the quark condensate is understood to be in the MS scheme. As 
shown in [21, 22], 788'' = —27^^ where 7^^ is the one-loop quark mass anomalous 
dimension^. This cancels exactly the scale /jq dependence in (3.7) to order q;s[21, 22]. 
The disconnected contribution to Im is thus 



2 .2r6T_ r-iFact _ q /nl^^ln\2^,, ^ l5'8(/fc)| 



2 



-^e^Fo^im cr* ^ ^o\m^^^c) ^^^^^ (4.15) 



but now with 



Ml 



Im CsilJiR) 



l + a(//c) (7S^ln^ + Ar88 



+ Im CMaiiin) (^^ In ^ + Argy) + 0{a^) . (4.16) 

Here one can see that the factorizable contribution is not well defined. It is due 
to the mixing of Q7 and Qs and is reflected here in the \vl(MxI ^J^r)- This Mx 
dependence cancels with the non- factorizable contribution of Qi in (4.9). Notice 
that the contribution of both terms, Im Cg and Im Cy, to Im are of the same 
order in 1/A^c. It is then necessary to add the non-factorizable term to have Im 
well defined. Since Im Gg is a physical quantity, factorization is not well defined for 
Qs. This was also shown to be the case for in [23]. Of course, the leading term of 
the l/A'^c expansion is well defined but that approximation would miss a completely 
new topology, namely the non-factorizable contributions. 

^For the other operators this correspondence does not hold and even for it is only valid in 
certain schemes, including ours. 

^See App. A for the explicit expressions. It can be seen there that no such relation holds for 

(2) 
788 . 
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4.2.2 The Connected Contribution 

From the leading high energy behaviour, the scalar-pseudo-scalar spectral functions 
satisfy in the chiral limit [24, 25] 

r dt - \lm uPsit) - Im U%{t)\ = = H dt- \lm ng^(t) - Im ngj,(t)" 

(4.17) 

which are analogous to Weinberg Sum Rules. Therefore the connected part of 
^^ss+pp(Q'^) satisfies an unsubtracted dispersion relation in the chiral hmit, 

nfs?;r\Q') = fit i 1^"^ . (4.18) 

Also in the chiral hmit, the scalar and pseudo-scalar (0 — 3) combinations satisfy 
other Weinberg-hke Sum Rules as shown in [26] for the scalar''' and in [27] for the 
pseudo-scalar, 

/ dt -Im ngj-') (i ) = = / dt -Im ngp^) {t). (4.19) 

Jo TT Jo TT 

We also know that the spectral functions Im Tiss{pp){Q'^) depend on scale due 
to the non-conservation of the quark densities. 

f^c^ Im ng(pp)(t) = 27^(/.c) Im ng(pp^(t) . (4.20) 

This scale dependence is analogous to the one of the disconnected part (4.15) and 
cancels the /ic dependence in Igsi/J'c)]'^ also for the connected part. 
We now proceed as for Qj and split the integral in (4.13) at /x^. 

4.2.2.a Qs'^""'' Long-Distance 

We perform simply the integral and obtain 

4.2.2.b Qs"'"" Short-distance 

Using the unsubtracted dispersion relation in (4.18), ^^ss+pp"^ iQ"^) ^ the chiral limit 
behaves at large in QCD as 



(Q^) -EE ^'^l^r^ (0|OJU)(0)|0>(-) (4.22) 

n=0 i=l H ^ ' 



''In [26] it was the alternative form of Eq. (3.12) which was used. 



10 



where 02^^_^_^^{Q) are dimension 2(n + 3) gauge invariant operators. 

Oi'\0) = O«(0) ; of (0) = of (0) . (4.23) 

Using the information on the mixing of Qi and Qg in (3.4), the scale dependence 
(4.20), it is easy to obtain the leading power behavior in (4.22) (see Appendix B) 

457r2 



2 

21l7r2 



a{uf + 0{a^); 
a{iyf + 0{a^). 



(4.24) 



Again the values of the coefficients of the power corrections in (4.22) can be 
calculated using global duality FESR, 



E E(-l)"^^'(0|^^(.n+3)(0)|0)(5o) 



m=0 i=l 



27ri Jcs, 



ds 



;l+m— n 



Mn+1 = 



Jo TT 



with n > 0, and sq the threshold for local duality for this two-point function. Again 
only terms with n — m survive at 0{as) and one gets 



167r2 V 
1 



dg" 



167r2 



dt 



Ml 



In 



t \ 1 



1^ 



4.3 Sum 



We now add all the contributions of Eqs. (4.4), (4.9), (4.15), (4.21) and (4.26) to 
obtain the full result. Notice in particular that all contributions contain the correct 
logarithms of Mx to cancel that dependence in Eqs. (3.6) and (3.7). 

The integrals over the spectral functions in the respective long and short-distance 
regime can in both cases be combined to give a simple ln(t//x^). 

Therefore, when summing everything to 0{as) and a// orders in we obtain 



-^e'Fo^Im Ge^ llu, C^{^lR) 



1 + ai^xc) I 777 In + Ar77 



-Mm C8{iiR)a{iic)^r7s 
+ I Im Cs(m) 



/I) 
9 

167r2 



Alr{iJ') 



1 + a{nc) I 788^ ln — + Args 



Im C7{[ij^a{[ic) 1^787 ^ + A'^s? j | 
X (3(0|gg|0)2(/.c) + j^^Sp(/^,/^c)) 



X 



(4.27) 



11 



where 




(4.28) 



To obtain this result we have used the local duality relations 




(4.29) 



The expression in (4.27) is exact in the chiral limit, at NLO in as, all orders in 
1/Nf. and without electromagnetic corrections. It doesn't depend on any scale nor 
scheme at that order analytically. The dependence on Mx also nicely cancels out. 
The lie dependence cancels against the fic dependence of the densities. Notice that 
in this final result we have taken into account the contribution of all higher order 
operators. 

This result differs from the ones in [3, 4, 5, 6] in the finite terms Arss and Aryg. 
They are necessary to cancel the scheme dependence. In addition, [4, 5] only take 
into account the dimension eight corrections and [6] uses a hadronic large Nc Ansatz 
to estimate them. 

As noticed in [12], the connected scalar-pseudo-scalar two-point function is ex- 
actly zero in U(3) symmetry, i.e. is l/N^ suppressed. We used this fact to disregard 
this contribution there. We will check later the quality of this approximation from a 
phenomenological analysis of its value. 

5. Bag Parameters 

We now re-express our main result (4.27) in terms of the usual definition of the bag 
parameters 



- e^Fo^Im Ge = {0\qq\Of if^c) [Im C^ifXRjBr^ifXc, m) + 3Im CsiliR) Bs^{fxc, f^n)] 
o 

= -6Im CM (0|O«|0)^(/x«) +Im Cs(iir) (0|Of |0)^(/x^) , (5.1) 



where the subscript x means in the chiral limit. 

This definition coincides with the one in [12] and gives 




12 



+ 3 a{fj,c) Arsr 



1 + 



1 



1 + 



« In ^ 



487r2 {0\qq\Oy{fic) 
I + Ar88 a{nc) 



X 



X 



1 + 



+ 



487r2 mqWiy^c) 
a{iic) . 3 



Aspil^R-, He) 



(5.2) 



The finite terms that appear in the matching between the X-boson effective the- 
ory with a cut-off and the Standard Model regularized with the NDR scheme were 
calculated in [12], 



Ar 



NDR 
77 



Ar 



NDR 
78 



NDR 
87 



1 

'8 



Ar 



NDR 



8 8Nr 



(5.3) 



The finite terms to pass from NDR to HV in the same basis and evanescent operators 
we use can be found in [10]. In the HV scheme of [9, 17] ^ these finite terms are 



A„HV-NDR 
^'77 



2N^ 



HV-NDR 
78 



1 Ar-^ 
) ^' 87 



HV-NDR 



3 
2 



Ar 



HV-NDR 



Nr. 



2Nr 



(5.4) 

The results for the scheme dependent terms Aryy and Ar87 [12] agree with those in 
[4, 6]. 

The -By and Bg bag parameters are independent of but depend on /ir and 
HC: and these dependences only cancel in the physical value of Im Ge- The jic 
dependence is artificial and a consequence of the normalization of the bag parameters 
to the quark condensate. 

At NLO in l/N^ we get 



B?^^^{^^c,^^R) 



1 + 



a{l^R) 
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24 ; 167r2 (0|gg|0)2(/xc) 
1 1 



Alr{iJ'r) 



1 + 



487r2 {0\qq\Oni^c) 



■Aspii^R, l^c) 



1 + 



1 

12 



54 In 



1 + 



1 



X 



9 



487r2 (0|gg|0)2(/xc) 



+ 23j a{jic) 



-ALpifil 



(5.5) 



647r^(0|gg|0)2(/xc)' 

*Thcrc is a finite rcnormalization from the HV scheme of [10] to the HV scheme of [9. 17], If one 
uses the Wilson coefficients in the HV scheme including the Cp terms from the rcnormalization of 
the axial current as [10], one has to add —Cp to the diagonal terms Aru in (5.4) and —(3iCf to 
the diagonal terms in the two-loop anomalous dimensions in (A. 7). 
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and in the HV scheme [9, 17] 



9 



33 , , 



1 + 



167r2 (0|gg|0)2(/xc) 
1 1 



X 



1 + 
1 + 



487r2 (0|gg|0)2(/ic) 
1 1 



X 



+ 



487r2 (0|gg|0)2(/xc) 



647r2 (0|gg|0)2(/xc) 



(5.6) 



We find an exact result for these B-parameters in QCD in the chiral hmit including 
the effects of higher dimensional operators to all orders. The scheme dependence 
is also fully taken into account. To our knowledge this is the first time these fully 
model independent expressions bag parameters are presented. 

6. The UlniQ^) Two-Point Function and Integrals over It 

There are very good data for Hj^ji [18, 19] in the time-like region below the tau lepton 
mass. They have been extensively used previously [4, 5, 28, 29], see the talks [20] for 
recent reviews. We consider it a good approximation to take this data as the chiral 
limit data. Nevertheless, one can estimate the effect of the chiral corrections with 
Cauchy's integrals around a circle of radius 4m^ of the type 



/ d5 5«in-(5)nL(-5) 



(6.1) 



with n > 0, m = 0, 1. For all the integrals we use, we have checked that these 
contributions are neghgible using the CHPT expressions for li^ji{Q^) at one-loop[24]. 
The discussion below is focused on the ALEPH data but we present the OPAL results 
as well. 

Wc reanalyze here the first and the second Weinberg Sum Rules (WSRs) [15], 
which are properties of QCD in the chiral limit 



(6.2) 



/ dt - Im Ul^it) = / dt - Im Ul^it) + Oia') = ; 

Jo TV Jo TV 

/ dU - Im n^o(t) = / dtt-lmUlJt) + 0(a'^) = 

Jo TT Jo TT 



where we used the perturbative QCD result for the imaginary part at energies larger 
than So, i-e. we assumed local duahty above Sq- These two WSRs determine the 
threshold of perturbative QCD sq- We used the experimental value for the pion 
decay constant = (92.4 ± 0.4) MeV. 
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Figure 1: The first and second Weinberg sum rule as a function of the upper integration 
variable s. The central curve corresponds to the central values of [18] while the upper and 
lower curve are the one sigma errors calculated as described in the text. 



These two sum rules are plotted in Fig. 1 for the central data values and the 
one sigma errors. These latter are calculated by generating a distribution of spectral 
functions distributed according to the covariance matrix of [18]. We then take the 
one sigma error to be the value where 68% of the distributions fall within. All errors 
in the numbers of this section and in the plots shown are calculated in this way. 

At this point, we would like to discuss where local duality sets in: sq. As we 
can see from Fig. 1 for s < there are two points where (6.2) are satisfied, the 
first one around 1.5 GeV^ and the second around 2.5 GeV^. Of course, this does not 
mean that local duality is already settled at these points as the oscillations show. 
One can expect however that the violations of local duality are small at these points. 
It is also obvious that local duality will be better when the value of Sq is larger. The 
procedure to determine the value of Sq is repeated for each of the spectral functions 
generated before and we use consistently a spectral function together with its value 
of the onset of local duality. 

There are several points worth making. Though for every distribution the first 
duality point in the 1st WSR is very near the corresponding one of the 2nd WSR, 
they differ by more than their error. Numerically, when used in other sum-rules they 
produce results outside the naive error. The second duality point, sq ~ 2.5 GeV^ 
yields more stable results. There is no a priori reason for the value of Sq to be exactly 
the same for different sum rules. 

Though the change from the 1st WSR to the second is small, and even smaller 
if one looks at negative moments, when one uses large positive moments (the ones 
we need here), the deviations are quite sizable as we will show. This is because pos- 
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itive large moments weigh more the higher energy region and the negative moments 
essentially use only information of the low energy region. 

Probably in the second duahty point, local duahty has not been reached either 
but certainly we should be closer to the asymptotic regime. We therefore take the 
highest global duality point available, the solution of Eq. (6.2), around 2.5 GeV^. 
Fortunately, for the physical matrix elements, the additional log(t//x^) in the inte- 
grand reduces the contribution of the data points near the real axis for t around 
/x^. This makes these sum rules much more reliable than the single moments used 



The second, and highest value with good data, value of Sq where the WSRs are 
satisfied runs roughly between 2.2 GeV^ and 3.0 GeV^. But not all of these values 
are equally probable. If we look at the distribution of the Sq values, there is a clear 
peak situated around the value calculated with the central data points but there are 
tails towards higher sq. The widths of the peak are essentially the same as the errors 
we quote. The sq where the second WSR are mainly in the area 



So = (2.53l°:J^) GeV^ (ALEPH) , so = (2.49^°;}^) GeV^ (OPAL) . (6.3) 



So = (2.56l°:J^) GeV^ (ALEPH) , sq = (2.53^°;}^) GeV^ (OPAL) . (6.4) 



These errors have been obtained as explained above. In the analysis below we use all 
experimental distributions with their associated value of Sq and not only those with 
So in the intervals above. 

The OPE of the n^^((5^) was studied using the same data[18] in [28]. They 
obtained a quite precise determination of the dimension six and eight higher dimen- 
sional operators from a fit to different moments of the energy distribution. This 
procedure has in principle smaller errors since one can use the tau decay kinematic 
factors which suppresses the data near the real axis but has a different local duality 
error. They use as upper limit of the hadronic moments, we agree with [29] that 
one should use the sq where there is global duality with QCD to eliminate possible 
effects of the lack of local duality at M^. Another comment is that as noticed in [6] 
the a'g corrections used in [28] are in a different scheme [7]. These corrections in the 
scheme used in [17] are presented in the appendices. 

We can determine the following higher dimensional operator contributions (4.5) 



in [4, 5]. 



and where the first WSR is satisfied in 




m=0 i=l 



oo 



;(m+3) 



(0)|0)W(-ir5^£ 



ds 



a 



2(m+3)(-^0' 
gl+m 




16 



E E(o|o?(U3)(o)|o)(3o)(-i) 



1 

2Td 



ds 



^2(L+3)("^0, 



(6.5) 



m=0 i=l 

In Fig. 2 we have plotted the value of M2 and M3 as a function of sq used in the 
integration, together with the one sigma error band. It is immediately obvious that 
the main uncertainty is the choice of sq to be used. This uncertainty is increasingly 
important with the increase of the moment. 
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Figure 2: The second and third moment as a function of the upper limit of integration 
So and the one sigma variation. 

Using ALEPH data on V — A spectral functions we get for the dimension six 
and eight FESR using the value for sq where the second WSR is satisfied 



M2 
Ms 



-(1.7+};^) ■ 10-3 GeV^ 
(7.2i|2).io-3GeV^ 



(6.6) 



The error bars are obtained by taking 68% of the generated distributions within this 
value, only including those where the WSR can be satisfied. The error is smaller 
than one would judge from Fig. 2 since the value of M2 and M3 at the value of sq 
where the spectral function satisfies a WSR is much more stable than the variation 
at a fixed value of sq. 

Using the OPAL data we get. 



-(2.0+};°) ■ 10-3 GeV^ 



M2 : 

M3 = (5.2^^:^) ■ 10-^ GeV« . 

Eq. (6.6) and (6.7) should be compared to the results in [18, 19, 28] 

M2 = -(3.2 ±0.9) ■ 10-3 GeV*^ 
M3 = -(4.4 ± 1.2) ■ 10-3 GeV^ . 



(6.7) 



(6.8) 
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The result for M2 is compatible within errors but M3 differs even in sign. Our error 
bars take into account the variation of sq but our result for M3 at the second duality 
point is always positive. 

This indicates a potential problem in the determination of M3 and higher mo- 
ments (and of smaller importance in M2). As said before violations of local duality 
can be sizeable for higher moments like M3 even at t ~ used as upper limit of 
the moment integrals. It would very helpful to do the same type of fit analysis done 
in [28] but using the duality point Sq. Our conclusion is that moments like M2, M3 
and higher are unfortunately unreliable unless one has data at higher energies. 

The integrals which are needed for Eq. (4.27) can be evaluated from the ALEPH 
data in the same way. We need 

AUfiR) ^ Tdtt^ln (-^] -Im ni^(t) = (4.7l°i) " 10"^ GeV^ 

A'^T'i^^R) = -f ° dtt' In (1 + ^) ^ Im nL(t) = (3.7t°i) ■ 10-=^ GeV^ ; 

A^lf^'M ^ r dt e In f 1 + 4) - Im nL(t) = (l.Otrr) " 10"' GeV^ (6.9) 



at = 2 GeV and using for each 
distribution its second duality point 
Sq. Notice the much smaller error of 
Alr and A^"^^^ when compared with 
M2 and M3. These values are all taken 
at the second duality point Sq where 
the second WSR is satisfied. We plot 
Alr SiS Si function of sq in Fig. 3. The 
OPAL data give instead 



AiRim) 



4 Lower 
•^LR 



A Higher 
•^LR 



(4.4 
(3.8+°-^ 
(0.6 



+0.4N 
0.3) 



0.5/ 

-0.8^ 
-0.6/ 



• 10"^ GeV^ ; 

■ 10"^ GeV^ ; 

■ 10-=^ GeV^ 

(6.10) 




Figure 3: The integral over the spectral func- 
tion needed for Im Ge- 



As a test, we can also calculate 
the electromagnetic pion mass difference in the chiral limit [31], 



B 



LR 



rso f t \ 1 

/ dU In ^ -Im n^R,-, ^ 

-(5.2 ± 0.5) • IQ-^ GeV^ (ALEPH) 
-(5.2 ± 0.6) ■ 10"^ GeV^ (OPAL) ; 



(6.11) 
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where wc also used the value of Sq given by the 2nd WSR. Notice that Blr does not 
depend on hr due to the second WSR (6.2). The experimental number is 

S = -(5.15±0.90) . lO-GeV. (6.12) 

where we used Fq = (87 ±6) MeV as the chiral limit value of the pion decay constant 
and removed the QCD contributions [32]. 

For comparison we quote the central values using as sq the second duality point 
where the first WSR is satisfied 

So = (2.56l[J:J:5) GcV^ , ALR{2GeV) = (4.0tg:?) " 10"'^ GeV^ , 

M2 - -(O.lllli) • 10-'^GcV^ ^i|j'"^"(2GeV) = (2.2 ±2.1) ■ IQ-'^GeV^ (6.13) 
M3 = (111?) • 10-3 GeV« , ^f^^''''(2GeV) = (1.812J) • lO"' GeV^ 

for ALEPH and for OPAL 

So = (2.531°:}D GeV^ , ^Li^(2GeV) = (3.4lgi) • 10"=^ GeV^ , 

M2 = (O.ltli) ■ 10-=^GeV^ ^|;?j^^^(2GeV) = {1.7t\i) ■ 10-=^GeV^ (6.14) 
Ms = (lOli) • 10-' GeV' , ^f^^'^^"(2GeV) = (1.7l}J) • 10-^ GeV^ 

The errors are larger here. The value of sq where the first WSR is satisfied varies 
more and is somewhat larger than the sq where the second WSR is satisfied, this 
makes the last results more dependent on the spectral function at high t which have 
large errors. 

If one tried to see the results using the first duality point, where less duality with 
QCD is expected, we get that using the one from the 2nd WSR 

So = (1.47 ± 0.02) GeV^ , ^^(SGeV) = (3.3 ± 0.1) • 10"=^ GeV^ 

M2 = -(6.6 ± 0.2) ■ 10-3 GeV^ ^i^"'""(2GeV) = (5.9 ± 0.2) ■ lO'^ GeV^ (6.15) 
Ms = -(121^) • 10-3 GeV^ >lf^^'^'''(2GeV) = -(2.6 ± 0.1) GeV^ 

Notice that M2 is not compatible with (6.6) with the central values differing by more 
than twice the error. The moment M3 changes even sign with respect to the second 
duality point showing the problems of local duality violations for larger moments 
more dramatically. As argued before one should the largest value of sq to ensure 
better local duality. However, the physical relevant moment Alr is much more 
stable with sq. 

7. The Scalar-Pseudo-Scalar Two-Point Function U^ss+ppiQ^) 

In this section we discuss some of the knowledge of the spectral function Im II^°^_^pp(t) 
which governs the connected contribution to the matrix element of Qs- In the large 
Nc limit there is no difference between the singlet and triplet channel so the inte- 
gral in (7.1) is 1/A^c suppressed and its contribution to Im Ge is NNLO. But in the 
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scalar-pseudo-scalar sector, violations of the large Nc behaviour can be larger than 
in the vector-axial- vector channel. It is therefore interesting to determine the size of 
this contribution as well. 

After adding the short- distance part to the long-distance part, the relevant in- 
tegral is (4.27) 

\ , , — - r dU In ( ^7 I - Im U^°^f>pp(t) , (7.1) 
487r2 (0|gg|0)2(/ic) io ^ ss+ppK J , I J 

This is the contribution of the connected part relative to the disconnected one. sq is 
the scale where in this channel QCD duahty sets in. The scale /x is the cut-off scale. 
The dependence on this scale being NNLO in 1/Nc cannot match the present NLO 
order Wilson coefficients. 

We can use models like the ones in [26] to evaluate the scalar part of the integrals. 
We only use the model there using the KLM[33] analysis, since only it ratifies (4.19) 
at a reasonable value of sq- In Fig- 4 we plotted for that parameterization the sum 
rule and the relative correction from the scalar part to the disconnected contribution 
3(0|gg|0)^(/i/j) for /i = /i/j = 2 GeV. The value of the scalar part of Eq. (7.1) is 



0.012 I 1 , 1 1 1 1 1.5 




-0.01 I 1 ' 1 1 1 1 -1.5 



0.5 1 1.5 2 2.5 3 0.5 1 1 .5 2 2.5 3 

s [GeV^] s [GeV^] 

Figure 4: The sum rule (4.19) as a function of sq for the scalar part and the relative 
correction to the disconnected contribution for Qg using the same parametrization as a 
function of sq for ^u/j = /x = 2 GeV. Notice that the correction is small in the region where 
the sum rule is satisfied. 

about 0.18 at sq = (1.41 GeV)^. 

In the large A^^^ limit, a sensible alternative estimate is to use meson pole domi- 
nance. In the pseudo-scalar sector, the U(3)xU(3) symmetry is broken by the chiral 
anomaly splitting the singlet rji mass away from the zero mass for the Goldstone 
boson octet. 
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Three meson intermediate states are not studied enough to be included at this 
level, we include instead the first tt' resonance. This means including a massless 
Goldstone boson plus the first tt' resonance and the singlet rji. The pseudo-scalar 
sum rule in (4.19) requires the following relation between the octet and the singlet 
couplings to the pseudo-scalar current for sq — 2.0 GeV^, 



m 



(7.2) 



Phenomenologically F^,/F^ << 1 [34] and Fq ~ F^^. 

We can introduce a scalar meson octet 5*8 and a singlet 5"! using the methods 
of [35]. The coupling constant for the octet can be denoted by Cm and has been 
estimated in [35, 36] to be about (43 ± 14) MeV. In fact, the sum rule (4.17) is a 
property of QCD and relates in this approximation Cm to Fq 



F' + Fi + 



8 



El 

8 



(7.3) 



which numerically agrees quite well with the phenomenological estimate. 

The scalar sum rule in (4.19) requires the singlet and the octet components to 
have the same coupling leading to a relative correction from the scalar integral to 
the disconnected contribution of 



F2 



127r2 F^ 



Si 



- Ml In 



M; 



(7.4) 



using both sum rules and the lowest meson dominance approximation. 

The contribution from the pseudo-scalar connected two-point function relative 
the disconnected contribution can then be evaluated to 



F^ 



127r2 F(f 



< In 



M 



(7.5) 



The contribution of the tt' is negligible. 

As said before the scale fi is free and cannot be at present matched with OPE 
QCD since it is a NNLO order in 1/A^c effect. The scale independence is reached 
when the sum rule 



Mi^ Ml- Ml- Ml 



p2 







(7.6) 



which is 0{N^a'^) (4.24) is fulfilled. This sum rule is very well satisfied in the linear 
(7 model, see e.g. [26]. 

The masses M^^ ~ 0.86 GeV (chiral limit value) M^/ ~ 1.3 GeV are known. 
The masses of the singlet and octet of scalars are not so well known. Using Ms^ — 
Ma ~ 0.5 GeV and Ms^ — MaQ(98o) ~ 0.98 GeV the correction to the disconnected 



21 



contribution is almost independent of the scale /i, neglecting the vr' and is independent 
of II for F^,/Fq = 0.017. The relative scalar contribution is about 0.40 and the 
relative pseudoscalar is —0.71 with a total relative contribution of about —0.31. The 
scalar contribution is within errors compatible with the earlier estimate. 

Prom the discussion here it can be seen that there is a possibly sizable correction 
but we expect it to be smaller than —40%. 

Of course, the scale dependence left in (7.1) is unsatisfactory in principle but 
small since the sum rule (7.6) is quite well satisfied. 

8. Numerical Results for the Matrix-Elements and Bag Pa- 
rameters 

The vacuum expectation value in the chiral limit of Qr itself is related directly to 
By^. This allows us to obtain^ 

(0|O«(0)|0)^^«(/.«) = - (l + la(/.«)) ^.AUf^n) 

(8.1) 



+ ^ aiim) 



mm\0)'{l^R) + Y^Asp{l^R,l^c) 



(0|O«(0)|0)^^(/.^) = - (l - lia(/x^)) ^ A«(//«) 



■ (8.2) 



For the numerics, we use the value of the condensate obtained in the MS scheme 
in [37], 

(0|gg|0)(2GeV) = -(0.018 ± 0.004) GeV^ , (8.3) 

the numerical results of Eq. (6.9), 

a(2 GeV) = 0.102 (8.4) 

and neglect, in first approximation, the integral over the scalar-pseudo-scalar two- 
point function. 

The weighted average of the first and second WSR results for ALni'^GeV) from 
ALEPH data is 

^t™(2GeV) = (4.5 ± 0.5) • lO^^ GeV^ (8.5) 
and from OPAL data 

-4°|^^(2GeV) = (4.2 ± 0.4) • 10"^ GeV^ . (8.6) 



^The analytical formulas are in agreement with [4, 5, 6] for the scheme dependent terms in Qj 
matrix elements but not for the Qs ones in [4, 5] and they were not included in [6] . 
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Though the systematic errors aver very correlated, sicne the central values are very 
similar we take the simple average of both results as our result 



and obtain 



and 



ALR{2GeV) = (4.35 ± 0.50) • lO'^ GeV^ 

(0|O^')(0)|0)^^^(2GeV) = -(4.0 ± 0.5) • 10-^ GeV^ 
= (-(4.2 ± 0.5) + (0.2 ± 0.1)) ■ 10"^ GeV*' 
= ((-3.3 ± 0.5) + (-0.9 ± 0.8) + (0.2 ± 0.1)) • 10"^ GeV^ 

(0|O^^^(0)|0)^^(2GeV) = -(6.2 ± 1.0) • 10"^ GeV^ 
= ((-3.9 ± 0.5) - (2.3 ± 0.9)) ■ 10"^ GeV^ 
= ((-3.1 ± 0.5) + (-0.8 ± 0.7) - (2.3 ± 0.9)) ■ 10"^ GeV^ 



(8.7) 



(8.8) 



.9) 



where we quote, namely, the total result, the integral and the vacuum expectation 
value separately and in the last case also the long and short-distance part of the 
integral separately. 

The short-distance part of the integral, the second term in the above, is the 
contribution of all higher dimensional operators. We find that its contribution is 
between a few % up to 35 % depending on the value of At fi — 2 GeV it is 
somewhat larger than the error on the integral cut-off at /i. 

Similarly, the matrix-element of Qg is directly related to and we obtain^° 



(O|0f (0)|0)^^^(2GeV) 
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1 + - a(2GeV) 



X 



X 



3(0|^?|0)^(2GeV) 



IGtt' 



;Asp{2GeV, 2GeV) 



27 
64^ 



a(2GeV) ALR{2GeV) . 



(8.10) 



(0|Of (0)|0)^^(2GeV) = 



l + -a(2GeV) 



X 



X 



3 (0|gg|0)2(2GeV) + ^^5p(2GeV, 2GeV) 



+ 



9 



647r2 



a(2GeV) ^LR(2GeV) 



(8.11) 



Using the same input as above we obtain 

(0|Of (0)|0)^-°^(2GeV) = (1.2 ± 0.5) • lO-'GeV' , 

(0|Of^(0)|0)^^(2GeV) = (1.3 ±0.6) • 10-=^GeV^ (8.12) 

^°We disagree in this case with the results in [4, 5, 6] because of the scheme dependent terms. 
These references also disagree with each other. 
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where the contribution of the integral over Im Hj^ji is at the 1% level and thus totally 
negligible. 

Another combination of these two matrix-elements can also be obtained from an 
integral over the ALEPH data[9, 17] by putting (4.7) and (4.5) in (6.5) including 
also the as correction of the appendix 

M2 = r dt - Im nL(t) = Y: («o, ^o) of {so) 
Jo n ^ 

>^^^(so) 

'6 

41 



13 



2{l + fa{so)) (010^(0)10)^ 



— a[so) 

+ (1 + f «(^o)) (0|Or(0)|0) 



NDR, 



So) 



An' 



a(so) 



2(l + -a{so)) (0|oW(0)|0)^^(.o) 



+ 
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aiso)) (O|0f (0)|0)^^(5o) 



= -47r^a(so) 



61 , , 



47 

1 + Y^«(so) 



167r2 



Alr(so] 



(8.13) 



The right hand-side is physical and we cheked that is independent of the scale Sq 
and scheme. We can therefore evaluate it at Sq = 4 GeV^. The contribution from 
(0|0^^^(0)|0)^^^(so) is numerically very small and we obtain 



Mo. 



-(2.0 ±0.9) • lO"'' GeV\ 



perfectly compatible within errors both with the result obtained from the data in Eq. 
(6.6) and with the result (6.8). This confirms our results on the size of the integral 
over Im Ilss+pp{t), which can therefore be considered negligible within the present 
accuracy of the disconnected contribution and M2. 

There is another sum rule which combines the two matrix elements, 



Ml 



^°dtt- Im nss+pp{t) = -E Cf {so, So) Of {so) 
^ i=i 



TT 



^-:-a{soy 211(0|Of (0)|0)(5o) + 90(0|0^^^(0)|0)(5o) +0{a'). 



(8.15) 

For the calculation of the coefficients see Appendix B. This sum rule is much less 
accurate than M2 since the leading terms are a| and the value of Mi is not known 
directly either. Therefore we don't use it. 

^^We thank Vincenzo Cirigliano, John Donoghue, Gene Golowich, Marc Knecht, Kim Maltman, 
Santi Peris, and Eduardo de Rafael for pointing out an error in the matching coefficients in the 
previous version of our paper. Our result agrees with the result found in [39] 
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The numerical estimates of the disconnected part, Asp, given above change these 
numbers somewhat but within the errors quoted. 

These results can also be expressed in terms of the bag parameters: 

S^^^(2GeV) = 0.75 ± 0.20 ; B^^{2GeV) = 1.15 ± 0.30 
s£^^(2GeV) = 1.2 ± 0.3 ; B^^(2GeV) = 1.3 ± 0.4 . (8.16) 

We can also express it in terms of Im Ge'- 

Im Ge = Im t(-2.1 ± 0.9) 10"^ GeV^ (8.17) 

which is quite compatible with the estimate in [12]. 



9. Comparison with earlier results 

To compare with other results in the literature we propose to use the VEVs (0|Og''|0) 
and (0|Og^''|0). The reason is that these quantities are what [4, 5, 6] and we directly 
compute. The matrix elements oi K ^ tttt through and Qs, in the chiral limit/^ 
can be expressed as follows [12] 



2^6 



{o\ol\o)M. 



x/6 (0|Of |0)^(/.«) 



3 Fi 



(9.1) 



The lattice results [38] are from computed K ^ tt matrix elements and use the 
physical values of fx and /j^ to convert into K tttt. Since we and [4, 5, 6] compute 
in the chiral limit, this amounts to a large spurious factor Jx/^/Fq ~ 1.6 of difference 
when comparing K tttt matrix elements or fKfn/FQ ~ 1.4 when comparing K ^ tt 
matrix elements. Usually this factor is not taken into account. Moreover each group 
uses different conventions, either the chiral limit value of and fx or their physical 
value. We give the lattice results for X ^ tt rescahng with the factor above. 

We agree with the old results by [4, 8] including the higher order operators. In 
fact, they also make an estimate of their contribution which agrees with our full 
calculation. We also agree reasonably well with [5]. We agree borderline with their 
new results[39] within errors, though their central value is almost twice ours for Qs- 



12 



See [12] for the definition of M2[Q7] and M2[Qs\- 
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Reference 


-105(0|O^'^|0)^^^GeV^^ 


103(0|O^'^|0)^^«GeV-^ 


Bt = Bg = 1 


5.4 ±2.2 


1.0 ±0.4 


This work (SS+PP=0) 


4.0 ±0.5 


1.2 ±0.5 


Knecht et al. [6] 


1.9 ±0.6 


3.5 ± 1.1 


Cirigliano et al. [39] 


2.7 ± 1.7 


2.2 ±0.7 


Donoghue et al.[4] 


4.3 ±0.9 


1.5 ±0.4 


Narison [5] 


3.5 ± 1.0 


1.5 ±0.3 


lattice [38] 


2.6 ±0.7 


0.74 ±0.15 


ENJL [12] 


4.3 ±0.5 


1.3 ±0.2 


Table 1: The values of the VEVs in the NDR scheme at /Xij = 2 GeV. 




Reference 


-105(0|O^'>)f^ GeV-^ 








Bj — Bs — 1 


5.4 ±2.2 


1.0 ±0.4 






This work (SS+PP=0) 


6.2 ± 1.0 


1.3 ±0.6 






Knecht et al. [6] 


11.0 ±2.0 


3.5 ±1.1 






Cirigilano et al. [39] 


8.2 ±0.9 


2.4 ±0.7 






lattice [38] 


4.3 ± 1.1 


0.8 ±0.2 






ENJL [12] 


7.1 ±0.9 


1.4 ±0.2 





Table 2: The values of the VEVs in the HV scheme at hr = 2 GeV. 



We do not agree with [6] even within errors. But if the O(a^) correction in 
(8.13) is taken into account, their result for Oq'^ goes to (2.7 ± 0.9) 10~^ and we are 
borderline also within errors though the central value is more than twice ours. 

We find a systematic a factor around 1.5 to 1.8 compared lattice results [38] for 
both matrix elements. They are compatible when we take the errors on both the 
lattice and our results into account and the fact that the lattice numbers are not in 
the chiral limit. The corrections for the latter we included only partially with the 
physical values of fx and /„-. 

We have not quoted the results from the CHPT large Nc approach[40] and the 
chiral quark model [41] because they are away from the chiral limit. 

Our results are very compatible with our earlier work. In [12] we used the 
ENJL model and a very low value of fi to estimate the same matrix-elements. The 
underlying reason for the agreement is that the operators and Qg mix quite 
strongly and the values of the matrix-element of Qs at the low scale fi ^ 0.8 GeV 
dominate the values of the matrix-elements of Q7 and Qs at the higher scale /i — 
2 GeV. The matrix element of Qs obtained here is the same as the one we used in 
[12], since the effect of Asp is estimated to be moderate here. 

9.1 Role of Higher Dimensional Operators 

We clarified here the role of the higher than six dimensional operators, an issue 
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raised in [8]. In our scheme they remove the /x-dependence which is not covered by 
the renormahzation group. 

The effect of higher dimension operators in our approach is to add A^'^^'^'^{^) to 
the low energy contribution A]^"^"^^ (fj,) , these are defined in Eq. (6.9), 

AUi^) = A'^f- + A^,T' = r A-^- nL(t) (9.2) 



where /i is an Euchdean cut-off. It is clear than the contribution of higher than 
dimension six operators is less important only for values of /i^ larger than sq, where 
Im n|^^ vanishes because of local duality. In Figure 5 we plot the two separate 
contributions and the sum as a function of /x. 

From the figure we can see for /i 
larger than 2 GeV the contribution 
of all higher dimensional operators is 
less than 25 %. We agree with [8] that 
for the matrix elements that involve 
integrals of Im one has to go to 
such values of /i to disregard the con- 
tribution of higher dimensional oper- 
ators. The contribution we find is some- 
what smaller than in [8] since we in- 
clude the effect of all higher order op- 
erators, not just dimension eight. 

The high value of is set by the 
threshold of perturbative QCD sq which 
depends very much on the spectral func- 
tion and on the integrand behaviour. 
In fact, from [42] one can see that 
relevant spectral function for the 27- 
plet coupling reaches the perturbative 
QCD behaviour very soon, from 0.7 GeV to 1 GeV. The OPE matched impressively 
well with the hadronic ansatz at such low values with just dimension six opera- 
tors. Therefore though higher dimensional operators appear one can expect smaller 
contributions in cases like (j27 and Re Gs- 

The matrix-elements studied in this paper might be special in the sense that 
they follow from integrals over spectral functions which have no contributions at 
short-distances from the unit operator or the dimension four operators. As the good 
matching at low scales in the example in [42] shows, the other quantities which have 
these contributions might have much smaller higher dimension effects. 




Figure 5: The separate contributions to Alr 
and the sum. "Higher" labels the effect of 
the higher than six dimensional operators in 
the short-distance contribution and "Lower" the 
long-distance part. 
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9.2 The Large Limit 

In the large Nc Umit, the spectral functions have only poles: 



- Im nL(t) = -F^S{t) - flMlSit -MD + Y: f^M^S{t - M^y) (9.3) 



TT 



V 



and 



-Im n 



TT 



(0-3) 
SS+PP 



it) = 4- 



8 E - Ml) - 8 E d,5s5(^ - Ml) 

Si Ss 



E K^it - Ml) + F^5{t) + E F^,S{t - 



(9.4) 



at a// values of t. The Weinberg-like Sum Rules, Eqs. (6. 2), (4. 17), (4. 19), assuming 
local duality holds above Sq, impose 



E flMl 



E fvMi 

V 

E fvM^- E /lAfl = o. 



A 

^1- 



f2. 

-^0 > 



(9.5) 



<S0 



i^o^+ E E K =8 E 



^1 



58 



8 E 

Si 



(9.6) 



We will use the expression (9.3) to get one of the two relevant integrals defined 
in (4.28) and the EM pion mass difference in (6.11) 



M2 <so 

^L«W = 2 E fvM^ln 

V 

M2 <so 



'Ml 



V 



Ml<so 



Ma' 



-2 E flM'^ln 



(9.7) 



We can also calculate the moments defined in (4.8) 



n+2 



M^<so M\<so 

E /^M?"-^^)- E flMT-"'^. 

V A 



(9.8) 



These sum rules and higher moments were studied in [29] with the MHA Ansatz^^. 
The results obtained using MHA are [29] 



MHA 



.3 - -(8.1±2.7) • lO-'GeV^ 
^™(2GeV) = (2.8 ± 1.0) • 10"' GeV^ , Bf/^ = -(5.9 ± 2.0) • 10-=^ GeVt9.9) 



^^Which in this case means that the spectral functions are saturated by the pion pole, the first 
axial- vector and the first rho vector resonances. 
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for any value of Sq larger than the 1st duality point, i.e. Sq > 1.5 GeV^ by construc- 
tion. 

These should be compared with our results (6.6) from the second duality point 
or (6.15) from the first duality point or the fit (6.8) [28]. The results using the MHA 
Ansatz are compatible with the ones using the first duality point within one sigma. 

If one uses the second duality point, the moment M2 agrees borderline within 
errors but their central value is more than twice our result (6.6). However the physical 
quantiy Alr{2GcV) agrees within errors with our second duality point though our 
central value is larger. The moment M3 is more problematic and we find that in the 
second duality point it changes sign. 

Also [28] obtained the values of the condensates of dimension six an eight which 
are proportional to M2 and M3 from a fit to different type of moments of the same 
data. The results are in Eq. (6.8). The moment M2 is compatible within errors with 
our results using the second duality point but not with the our results using the first 
duality point. In particular, the moment M3 is incompatible with both our results 
using the first and second duality points. This must be due to the duality violations 
being weighted differently in (6.8) as mentioned previously. 

We conclude from these comparisons that for the dimension eight moment M3, 
local duality at low values of Sq starts to be a problem. As we go up in the moments 
we need more and more accurate information at higher energies. This affects also to 
M2 . 

The case of more than one resonance in each channel was analyzed in full gen- 
erality in [43]. There one can find the large differences that the p' produces in the 
second, third moments, and higher moments with respect to the case of a single 
vector resonance in each channel. 

10. Conclusions 

In this work, we have calculated in a model independent way the matrix elements of 
the AS — 1 operators Q7 and Qs in the chiral hmit. We have done it to all orders 
in 1/A^c and NLO in ag- 

The scheme dependence has been taken into account exactly at NLO using the 
X boson method as proposed and used in [11, 12, 23]. In fact, these two operators 
are a submatrix of the ten by ten done in [12]. 

We would like to mention some issues sometimes mixed up in the literature. First, 
the X boson method has nothing to do with using or not large Nc- It can be used 
without the large Nc approximation as well, as shown again in this paper. Second, 
our method of treating the scheme dependence is consistent and we never mix up two 
different schemes, cut-off and MS schemes. We do an analytic matching between 
a cut-off regularization and dimensional regularization in a well-defined scheme at 
perturbative scales first. The finite parts arising in this matching appear in the 
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methods using dimensional regularization to the end as well as explained in the 
appendix. 

Wc obtain exact matching in an Euclidean- cut- off regularization and analytical 
cancellation exact of (all) infrared and UV scheme dependences. 

For the contribution of higher order operators discussed in [8] and [5] we clarify 
how to include all higher dimensional operators and exact scheme dependence at 
NLO in as of both the Q-j and Qs niatrix elements. As a result we find smaller 
corrections due to this effects as discussed in Section 9.1. In our approach the effect 
of the higher order operators is to remove the remaining dependence on the Euclidean 
cutoff ji beyond the RGE cvohition. The result of resumming all higher dimensional 
operators in the case of Q'j makes its prediction much less sensitive to the choice of 

So- 

As noticed in [6, 12], Asp is zero in the large limit and therefore is Zweig sup- 
pressed. We find no sizeable violation of the dimension six FESR using factorization 
for Qs- 

We find that the moment M2 is very sensitive to the spectral function arount 
2 GeV2. 

Our main analytical results are the expression for the matrix-elements (4.27), the 
bag parameters (5.5), (5.6) and the expansion coefficients of the spectral functions 
(A.24), (A. 25) and (B.12). The main numerical results are the VEVs (8. 8), (8. 9) and 
the bag parameters (8.16). These results are exact in the chiral hmit, so we have the 
A7 = 3/2 part of e' /e model independently at all orders in 1/A^c- In order to reach 
final values all effects which vanish in the chiral limit, as final state interactions, 
quark-mass effects, isospin violation and long-distance electromagnetic effects still 
need to be included. 
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A. Calculation of the Corrections of 0{a^) to the Dimension 
Six Contribution to nJ^(Q^) 

A.l Renormalization Group Analysis 

We have the two-point function 

nr«(?) ^1^1 d^ye^^-y{Q\T{L>^{y)R''{OnO) ^ {q^q'' - g^'q'M^iq') 

+ ?Vni^(g2) (A.l) 

The contribution of dimension six operators to ^J,r{Q^) (where Q"^ — —q^) can be 
written in D = 4 — 2e dimensions as 

Q' nL(Q') ^ E Q") <Oi> (^) (A.2) 



D=6 

i=l,2 



with 



< Oi > = (0|OW|0) = i < 0|(s7^d)i(d7,s)R|0 > 

< O2 > = (0|O^^^|0) = 3 < 0|(5d)L(ss)K|0 > (A.3) 



and 



C,{v,s) = a{v)Y.a{vfc\'\v,s) (A.4) 

fe=0 

where the dependence in v and s of Cf^ is only logarithmic. Everything here we 
define in the MS scheme. 

In absence of electromagnetic interactions the matrix elements (A.3) only mix 
between themselves. The renormalization group equations (RGE) they satisfy are 

= -l77{i^) <0i> (z/) + ^787('^) < O2 > (z/) 

J < O2 > (ly) ^ ^^^^ + 6778(i^) < Oi > (ly) (A.5) 

au 

With 7(i/) the QCD anomalous dimension matrix defined in (3.4). In the NDR 
scheme [9, 10, 17] for n/ = 3 flavours^^, 



7(^) = j:^^-^a{ur 

n=l 

' 2N. \ Nr. 



NDRi2) ^ \_ ( -13riV2 + 132iV, - 45 



213Ar3 - 72NI + lOSA^e \ 
203Ar4 - 60Ar3 _ 479jy2 ^ ;L32Are - 45 J ' 

(A.6) 



^^For these operators the Fierzed version and the QtQs version have the same anomalous di- 
mension matrix. 

^^We will use along this work n/ = 3 since this is the number of active flavours of the QCD 
effective theory where Qr and Qs appear. 
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In the HV scheme of [9, 17] 



16 



7 



HV{2) 



1 / -17N^ - 12N^ - 45 



llbN^ - 12N^ - 71N^ - 12Nr - 45 



(A.7) 



We also need the quark mass anomalous dimension in the MS scheme, 
where m is a quark mass. The first coefficient is scheme independent 



(A.8) 



F ■ 



(A.9) 



Notice that 738^ = ~'^'ym all orders in l/N^ [21, 22], this is the reason why in 
the chiral limit is very near to 1 [12]. The large Nc result absorbs all the one-loop 
scale dependence. This exact scale cancellation does not occur for Qq even at leading 
order in as- There is a remnant diagonal anomalous dimension at one- loop of order 
one in 1/A^c which is not taken into account by the large Nc matrix clement. There is 
therefore no reason to expect Bq around 1 as sometimes is claimed in the literature. 
7^'' is the same for both the NDR and HV schemes[44]. 



M5(2) 
im 



203iV^^ - - 9 



(A.IO) 



The relation 738^ — —2'^^ is not valid: 

AT. 
788 



NDR{2) ^ _2^M5{2) _^ ^ 



32Ar2 



89A^; - 24A^e + 18 



(A.11) 



The two-point function ^^lr^Q^) independent of the scale u in. D — A 



D=6 



0. 



(A.12) 



This is also true in D dimensions if 75 is anti-commuting like in the NDR scheme. 
The HV results are obtained from the NDR ones using the published results in [10]. 
In D = 4 — 2e (A.12) yields the general condition 



= «'(^) (/^(«)(^ + 1) - 2eA;)Cf V, Q') <0i> (i^) 



k=0 
+ V 



dv 



(A.13) 



^^I.e. without the /3i Cp terms from renormahzing the axial current in the diagonal coefficients 
[17]. 
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with 



da{v) 
dv 



a {I3{a) - 2e) 



(A.14) 



and /3(a) = Y^k=i Pk^ii^)'^ with first coefficient Pi — 1 — llNc/6 for Uf — 3. 
To order a{u)^, one gets 



dv 



= 



(A.15) 



so the Cf^ are constants. 
To order a{y) 



du 

(0) c?CfM!) 
dv 



0, 



(A.16) 



Integrating these two equations we obtain 



with 



(1) 



2e 
2e 



+ 



(1) 



2e 
2e 



(1) 



(1) 



(1) 



a 



(0) 



^(0) .,(1) 



(A.17) 



(A.18) 



which are valid in D = 4 — 2e. The coefficients Cf'\ Df \ and f/^'* depend on e. The 
anomalous dimensions j3i and do not depend on e in MS", and in MS schemes in 
a known fashion. 

A.2 Calculation of the Constants Cf^ and f/^^ 

The bare vacuum expectation value of < Oi > can be expressed as an integral as 
follows 



< Oi > 



bare 



i f d^q 
7,9 



(27r) 



D-l 



/|^(Q^nL(g^)). (A.19) 



The scheme used here to regularize this integral is the MS scheme with D = 4 — 2e, 

dQ\Q'Y-%Q^Tll^{Q^)) (A.20) 



' 327r2 r(2 - e) Jo 

Notice that < Oi >bai"e jg scale independent. The integral (A.20) diverges due to the 
high energy behaviour of n|^^((5^). It is enough then to use the large expansion 
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of Uj^j^lQ^) in D dimensions. This is a series in (l/Q^)" starting at n = 3 in the 
chiral hmit, Eq. (A.2). Each coefficient of this series is finite and can be written as 
a Wilson coefficient times the vacuum expectation value of some operator. We now 
put (A.2) and (A. 17) in (A. 20) and perform the integral to find the divergent part. 
For that we need the integral, 



{Q 



2U+e 



e 



-2e 



(A.21) 



We will set fi = u afterwards. 



The MS subtraction needed then gives the full dependence on u. 



dv 



+ 



167r2 
3 

T67r2 



a{v) { Cf + a{u) 



+ 



D 



(1) 



< Oi (z/) 



< O2 (u) . 



(A.22) 



Overlined quantities are in four dimensions and 



(1) 



hm - ( - D 



(1) 



(A.23) 



Comparing (A.22) and (A. 5) order by order in a and using (A. 18), we get up to the 
needed order in e 



CI 



(0) 



IGtt^ 



■ (1) , 

777 +P77e 



a 



(0) 



9 



■ (1) I 

787 + P87 e 



(A.24) 



D«,D«,G«andG(^) are then determined up to the pij from Eq. (A. 18). We also 
get 



(1) 



167r2 



(2) 
777 



i:D« - ; 
6 ^ 2 ^ ' 



Stt^ 



(2) 
787 



l-o(i) 



ID';' - ^G^^^ 



1 

-( 
2 



(A.25) 



9 '°' 6 
The constants pij we determine below. 

A. 3 The constants pij 

We now evaluate Eq. (A. 20) to 0{a) fully with its subtraction in dimensional reg- 
ularization using the same split in the integral at //^ as we used in the main text. 
The short-distance dimension six part is the only divergent part, now regulated by 
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dimensional regularization rather than the X-boson propagator as in the main text. 
The result is 



^ Alr{i^). (A.26) 



327r2 



Comparison with Eq. (8.1) and (8.2) allows to determine P77 and psr- The finite 
coeflicients there are basically the Avij that corrected for the dimensional regular- 
ization to the X-boson scheme. If one works fully in dimensional regularization, it 
is here that these finite parts surface. 
The result is 

^NDR _ 3 „NDR _ 3 _ 

rf/ = -j|;; P?/ = J. (A.27) 

The transition between both agrees with the results in [10]. 
Putting in (A. 6) and (A. 7) to obtain numerical values 

^F^(o) _ 1 . ^-n(i) - _^ . 

IGtt^ ' 2' 167r2^ 2' 

l,D?^l (A.28) 



87r2 ^ 2' 87r2 

For (A. 25), in the NDR case we get 



3 -^NDR{1) 13 9 -^NDR{1) 75 



167r2 ' 16 87r2 ^ 16 

and in the HV scheme 

_ 3 ^ 41 9 ^ _63 

167r2 ^ 16' 87r2 2 15 ^ ^ 

All the expressions above are for n/ = 3 flavours. 

B. Calculation of the Corrections of O(a^) to the Dimension 
Six Contribution to Il^gg^^pp^^{Q'^) 

B.l Renormalization Group Analysis 

The function we have to study here is 

nSJ'p(g) ^ ^ / d^ye^y^ {Q\T[{S + tPY'-'\y){S - zP)(°-^)(0)]|0) (B.l) 
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with the definitions appearing in Section 3. 

The contribution of dimension six to the connected part of ^^ss+ppiQ^) ^® 
written as (4.22) 



i=l,2 



(B.2) 



with 



k=0 



and the operators Oi and O2 were defined in (A. 3). 

From (4.20) and (4.18), we have now in D = 4 — 2e, 



(B.3) 



(B.4) 



Using this relation and the renormahzation group equations, we get 



(1) 



2e 



+ 



2\ 



(1) 



2e 
2e 



(1) 



(1) 



(B.5) 



with 



.(0) 



(1) _ j^y 



(0) 



Pi - 27 



(1) _ j^y 

m 788 



+ g<^l 787 ■ 



(B.6) 



In the next Section we determine the values of the constants Cf^^ and F^''' , which 
depend on e. 

B.2 Calculation of the Constants ^ and 

The connected part of Il^gg^^pp{Q'^) can be related to the bare vacuum expectation 
value of the connected part of < O2 > i^^) through the relation 



< >^oTn (^) 



d Q -pf (0-3) conn/ 



(27r) 



D '-'-SS+PP 



{^) = j ^^'sAT\Q') (B.7) 



In the MS scheme with D = 4 — 2e and with renormalized li^gg^pp^Q"^) 



< O2 > 



bare 



(B.8) 
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Proceeding analogously to the case of n^^(Q^) in Appendix A. 2 and using that 
there is now a non vanishing contribution coming from the anomalous dimensions of 

n(0— 3) conn i 
ss+pp > namely, 



(27r)^ du 



= 27^(z/) < O2 > 



bare 



(B.9) 



that we have to add to the one from the i/-dependence of the subtraction determined 
by the integration of Q'^ in (B.7). 

The scale dependence of the total < O2 > can be obtained by adding both, we 
get in = 4 - 2e 



d < O2 >^'^ jiy) 
dv 



87r2 



a{v) < 



^(0) 

Ci +a(z/) 



^(0) 

C2 + a(i^) ' 



G« -.(1) 



+ 



^(1) 



Again the barred quantities have to be taken at e = and 



G. 



(1) 



hm - ( 



(B.IO) 



(B.ll) 



Comparing this equation with (A. 5) order by order in a{v) one obtains. 



(7{°^ = 487rV8e; 



7S^ + 27i'^+P88e' 



(2) 1^(1) if;«. 



^(1) 
2 



'-G 



(1) 



1^(1) 



(B.12) 



2 " 2 

Using Eq. (B.6) everything can then be determined in terms of the pij. 
B.3 Calculation of the pij. 

We now evaluate also the finite part from Eq. (B.l) fully in dimensional regulariza- 
tion to 0{a) and obtain 



C;'+487rV8) <Oi (u) 



+ iC,' -8n%s] <02 {u) 



+ 3 < 0|gg|0 >^ (u) 



(B.13) 
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Comparison with Eq. (8.10) allows to determine pyg and pss- The finite coefficients 
there are basically the Avij that corrected for the dimensional regularization to the 
X-boson scheme. If one works fully in dimensional regularization, it is here that 
these finite parts surface for the Qs contribution. 
The results are 

NDR _ _^ TV L • n^DR _ 3 . 

P88 - , P78 - 2 ' 

= pf3^ = -l. (B.14) 

again agreeing with the transition between both from [10]. 
Putting numbers, we get 

^(0) ^(0) ^(1) ^(1) 

= C2 =D, =D^ =0 

48^^^^ -32' "8^^^^ ^^-^^^ 
which are scheme independent. All the expressions above are for n/ = 3 fiavours. 
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